Estimation of Scale and Hurst Parameters of 
Semi-Selfsimilar Processes 

S. Rezakhah* A. Philippe* N. Modaresi* 



Abstract 

The characteristic feature of scmi-sclfsimilar process is the invariancc of its finite 
dimensional distributions by certain dilation for specific scaling factor. Estimating the 
scale parameter A and the Hurst index of such processes is one of the fundamental 
problem in the literature. We present some iterative method for estimation of the scale 
and Hurst parameters which is addressed for semi-selfsimilar processes with stationary 
increments. This method is based on some fiexible sampling scheme and evaluating sam- 
ple variance of increments in each scale intervals [A"~^, A"), n G N. For such iterative 
method wc find the initial estimation for the scale parameter by evaluating cumulative 
sum of moving sample variances and also by evaluating sample variance of preceding 
and succeeding moving sample variance of increments. We also present a new efiicient 
method for estimation of Hurst parameter of selfsimilar processes. As an example we 
introduce simple fractional Brownian motion (sfBm) which is semi-selfsimilar with sta- 
tionary increments. We present some simulations and numerical evaluation to illustrate 
the results and to estimate the scale for sfBm as a semi-selfsimilar process. We also 
present another simulation and show the efficiency of our method in estimation of Hurst 
parameter by comparing its performance with some previous methods. 
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1 Introduction 



Self-similarity has been discovered, analyzed and exploited in many frameworks, such as nat- 
ural images [6] , fluctuations of stock market [7] and traffic modeling in broadband networks 
|19j . The most frequently identified properties of high resolution traffic measurements from 
a wide range of packet networks are long-range dependence and self-similarity [22]. The 
fractal behavior of some important processes, such as fractional Brownian motion (fBm) and 
fractional Gaussian noise (fGn) is described by a single parameter, called Hurst parameter. 
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Self-similar processes are stochastic processes that are invariant in distribution under 
suitable scaling of time and space. These processes have increments which enter naturally 
in the analysis of random phenomena (in time) exhibiting certain forms of long-range de- 
pendence [16]. Semi-selfsimilar or discrete scale invariant processes requires invariance by 
dilation for certain preferred scaling factors only. In practice detecting self-similar property 
and estimation the self-similarity index and also scale parameter of the semi-selfsimilar pro- 
cesses are the main objects in the study of such processes. There are many problems arising 
in estimation of the Hurst parameter H and determine those effects which can influence the 
results considerably. Estimation depends on several factors, e.g, the estimation technique, 
sample size, time scale, level shifts, correlation and data structure. Beran [7], Taqqu and 
Teverovsky [27], Bardet et.al. [5] give a good review of statistical aspects of parameter 
estimation methods for self-similar and long-memory processes. Among these the most well 
known are variance and covariance based methods (rescaled adjusted range R/S statistic 
and variogram including aggregated variance). The other methods are maximum likelihood 
based methods. Whittle estimator, the local Whittle [25], the local log-periodogram estima- 
tors, the global log-periodogram estimator and roughness-length method j30j . Traffic model 
based on the fBm contains three parameters: the mean rate, variance parameter and Hurst 
parameter. Coeurjolly [13] developed a class of consistent estimators of the parameters of a 
fBm based on the asymptotic behavior of the /c-th absolute moment of discrete variations of 
its sampled paths over a discrete grid of the interval [0, 1]. He derived explicit convergence 
rates for these types of estimators, valid through the whole range Q < H < 1. Consistent 
estimators of the fractal dimension of locally self-similar Gaussian processes based on con- 
vex combinations of sample quantiles of discrete variations of a sample path, almost sure 
convergence and the asymptotic normality for these estimators are derived |14j . 

Several such studies of variations uncovered a generalization of fBm to non-Gaussian 
processes known as the Rosenblatt process and other Hermite processes. The processes are 
of order q and with Hurst parameter H G (^,1) which are self-similar with stationary in- 
crements and exhibit long-range dependence. The variations of these processes are studied 
and a consistent estimator for the self-similarity parameter from discrete observations of a 
Hermite process is constructed [TT] . For identifying a locally self-similar Gaussian process a 
new approach is proposed based on the asymptotic behavior of convex rearrangement, i.e., 
sums of increasing ordered variations [21] . They stated their result concerning the construc- 
tion of the estimators of the local Holder index. For the estimation of the Hurst parameter 
six methods are used that can be classified as temporal, spectral and time-scale methods, 
respectively. The temporal methods selected are (1) Rescaled range analysis (Beran [7]), 
(2) Level of zero crossings (Coeurjolly [12]) and (3) Detrended fluctuation analysis (Hu et 
al. [IT]). From the group of spectral methods, the (4) log-periodogram method (Beran 
[7]) is used while the (5) Wavelet transform modulus maxima method (Arneodo et al. [3]) 
and (6) Abry-Veitch estimator (Veitch and Abry [28]) are chosen from the group of the 
time-scale methods. Recently, moving average method to estimate the Hurst exponent and 
the correlation properties of the time series are presented and is compared with the rescaled 
range method of Hang Seng Index data for some periods [29]. 

Of late there has been an increasing interest in describing and discussing the integer 
valued analogues of classical distributions, like stable, semi-stable, semi-selfdecomposable 
and geometrically infinitely divisible [2] , [9] and [10] . The idea of binomial thinning to arrive 
at the right definitions of non-negative integer valued (/o-valued) semi-selfsimilar processes 
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are introduced |26j . In this class, they characterized semi-selfsimilar Levy processes in terms 
of an /o-valued first order autoregressive series. 

The paper is organized as follows. In section 2, we present flexible sampling method 
and some basic notions related to the self-similar and semi-selfsimilar processes in discrete 
parameter space. We introduce an example of semi-selfsimilar processes, named, simple 
fractional Brownian motion (sfBm) in this section too. We present some simulations for 
sfBm with different scales and also simple Brownian motion with different Hurst indices 
as semi-selfsimilar processes in section 3. In Section 4 we consider an iterative estimation 
method for scale parameter of semi-selfsimilar processes with stationary increments based 
on sample variances of increments. For this iterative method we find the initial scale pa- 
rameter by two procedures. In the first, we use the cumulative sum of moving sample 
variances (MSV) and in the second, sum of preceding and succeeding sample of variances 
of increments. We also present a new estimation method for estimating the Hurst param- 
eter of selfsimilar processes based on sample variance of different increments and show the 
efficiency of such new method by comparing its performance with some previous methods 
via simulation. 



2 Method of flexible discrete sampling 

We consider certain flexible sampling scheme as sampling at equally spaced points in 
each scale interval. Following this method of sampling from a semi-selfsimilar process 
{X{t),t G M"'"} with scale A > 1, we decide to have some flx number of samples in each 
scale, say T. So we provide a discrete time semi-selfsimilar process X{-) with parameter 
space {A"-i/c^^, n G N, A: = 0, 1, • • • T - 1}. 

A process {X(t),t G R^} is said to be self-similar of index H > 0, ii for every A > 0, 

{X{\t)} ^ {A^X(t)} (2.1) 

where = is the equality of all finite-dimensional distributions. As an intuition, self-similarity 
refers to an invariance with respect to any dilation factor. However, this may be a too strong 
requirement for capturing in situations that scaling properties are only observed for some 
preferred dilation factors. This process is said to be semi-selfsimilar with index H and 
scaling factor Aq > if (2.1) holds just for A = Aq, A; G N. 

A process {X{k), A: G T} is called discrete time self-similar process with parameter space 
T, where T is any subset of countable distinct points of positive real numbers, if for any 

^ ^ {Xik,)} ^ (^f{x{h)}. (2.2) 

The process X{-) is called discrete time semi-selfsimilar with scale / > and parameter 
space T, if for any A;i,A;2 = Ih G f , (2.2) holds; see [H]. If the process {X{t),t G M+} 
is semi-selfsimilar with scale / > 1, then by the above method, one can find a by the 
equation I = , T G N, then by sampling of the process at points a^,k G Z, we have 
X{-) as a discrete time semi-selfsimilar process with parameter space T = {a'^,A: G Z}. 
Based on the definition of wide sense self-similar process presented in [23], a random pro- 
cess {X{k),k G T} is called discrete time self-similar in the wide sense with index H > 
and parameter space T, where T is any subset of distinct countable points of positive real 
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numbers, if for all k,ki E T and all c > 0, where ck, cki E T : 

(i) E[X^{k)] < oo, 

(n) E[X{ck)] = c"E[X{k)], 

{Hi) E[X{ck)X{cki)] = c'^^E[X{k)X{ki)]. 

If the above conditions hold for some fixed c = cq, then the process is called discrete 
time semi-selfsimilar in the wide sense with scale cq. Through this paper we are dealt with 
the wide sense self-similar and wide sense scale invariant processes, where for simplicity we 
omit the term "in the wide sense" hereafter. 

Using quasi Lamperti transform as CH,aX{t) = t^Y{\og^t) |21j . the counterpart of any 
self-similar process X[-) can be identified as discrete time stationary process 

y(n) = /:^;„X(n) = a-"^X(a"). 

It is clear by this relation that \i X{-) is a discrete time semi-selfsimilar process with scale 
/ = , T € N and parameter space T = {a^,k € Z}, then Y{-) is a discrete time 
periodically correlated process with period T and parameter space T = {n,n € Z}. 

2.1 Simple fractional Brownian motion 

As an example of a semi-selfsimilar process, we introduce a process X[t) with index H > 
and scale A > 1 which is called simple fractional Brownian motion (sfBm) and is defined as 

oo 

X{t) = A("-i)(^-^')/[,„-i,,„)(t)i?H'(t) (2.3) 

n=l 

where Bh'{-), /(•) are fractional Brownian motion with index H' which is selfsimilar process 
and indicator function respectively. For H' = ^, B}ji{t) is Brownian motion and X{t) is 
called simple Brownian motion (sBm) which has Markov property [21] • It is easy to find 
that X{t) is semi-selfsimilar with scale parameter A but not Markov. 

3 Simulation 

In this section we have generated and plotted simple Brownian motion (sBm) defined after 
equation (j2.3p with scale A = 1.2 and Hurst indices H = 0.3, 0.5, 0.7 at points A"A;(A — l)/T 
for T = 20 and n = 1, . . . M, and A: = 0, 1, • • • T - Iwhere M = 20. Note that the fractional 
Brownian motion with index H' is simulated using the circulant matrix embedding method 
(see [5 for a description). 

We have samples of M = 20 scale intervals [A",A"~^"'^) and in each scale interval we 
have T = 20 equally space samples. Simple Brownian motion is a semi-selfsimilar process 
with index H, and for H = 0.5 it will be Brownian motion, which has been plotted for 
comparison. By this method we show that how samples are effected by Hurst index. 

We simulate and sfBm as a semi-selfsimilar process with H = 0.9, and Hurst index of 
the fBm as H' = 0.7 for different scale parameter A = 2, 4, 8. 

Figure 1 (Left) is the plotted sfBm which shows how the process is effected by scale index. 
Both sBm and sfBm have semi-selfsimilar property with stationary increments. In Figure 2, 
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Fi gUrC 1: sfBm with Hurst indices H'=0.9 and H=0.7 and scales 2 [blue], 4 [green] ,6 [red] and iBm with Hurst 
indice H=0.7 [black]. (Left), sBm with scale 1.2 and Hurst indices 0.3 [blue], 0.5 [green], 0.7 [red] (Right), both with 
geometric sampling and equally spaced in each scale interval 
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Figure 2: Moving average Sample variances with equally space sampling of sfBm, Left figure Scale is 2 and right 
figure scale is 4. 

we show how our method for recognizing scale parameter works. We have simulated sfBm 
at equally space points with H = 0.9 and H' = 0.2 corresponding to equation (2.1) for 
scale A = 2 and scale A = 4 left and right correspondingly. We also evaluated and plotted 
corresponding moving sample variance (MSV) of some small number of observations, say 
b* = 10, which is defined as 

^. = ^E(nt.+.)-m)', (3.4) 
i=i 

where Y(ti) = ^ Sj=i ^(^i+j)) i = 0,1, ... ,n — b* , and n is the number of observations, 
in Figure 2. This figure shows how the sample variances can be clustered for scale intervals 
and have jump at starting point of each scale intervals. Then the mean of the ratio of the 
length of successive scale intervals provides the initial estimate of scale parameter for the 
iterative estimation method which is described in section 4.1. 
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4 Estimation procedure of parameters 



As semi-selfsimilar processes are characterized by certain scale parameter A and Hurst 
index H. For the estimation of scale, one need to consider samples of several scale intervals 
(A''~^,A'^], € N. For Estimation of the scale and Hurst index, we present a heuristic 
iteration method in this section which is proposed to estimate the scale parameter first and 
then the Hurst index. Our method can be applied in the situations where the process has 
stationary or independent increments. For such iteration method we need some starting 
points which has been described in sectionl. 

4.1 Initial choice of scale parameter 

The starting point of our estimation method can be evaluated by the following methods. 
Suppose some equally space sampling of the process with small enough space to have a 
large number of samples, say A^. Using (3.4), one could evaluate moving sample variances 
Vi for some small number, say 5*, of samples and plot Vi against i, for i = 0, 1, . . . , iV — 6*, 
see Figure 2. This plot provide some clusters for the sample variance where jumps occurs 
at the starting point of all scale intervals [A'^"^, A''"), A; G N. To cause these clusters to be 
distinguished well enough, we consider moving average of sample variance ViS defined by 



for some d* G N say d* = 20, and i = 1, 2, • • • n* where n* = n — b* — d*, b* defined by (3.4). 
By choosing appropriate values for bj and d*, and plotting consecutive samples of Wi's one 
could provide a clear distinguished clusters for samples of consecutive scale intervals. By 
detecting start points of these consecutive scale intervals, one could estimate an appropriate 
initial value for the scale parameter A. By calculating the cumulative sum ( Wj's ) defined 
by (|4.ip as Ui, . . . ,Un, where Ui = YTj=i ^j-> it is shown for example 1, by Figures 5(a) 
and 5(d), the plot of the sequence Ui provide some broken lines that broken points occurred 
at the border of the corresponding clusters, or scale intervals at Figures 4(c) and 4(f). So 
we can detect the starting points of scale intervals as these points and by these, the scale 
parameter of a semi-selfsimilar process. We consider following methods for such estimation. 

First Method: 

We can detect the starting points of successive scale intervals as change points on the mean 
of the sequence Wi. Let Tj, z = 1,2,3 to be the last three change points, or equivalently 
starting points of the last three scale intervals. So the initial choice of the scale parameter, 

say Ao, can be evaluated as — ^. 

T2 - n 

Many methods for finding single or multiple change points are available and implemented 
in the statistical softwares like R with the package "changepoint". Since the amplitude of 
the jump and the variance of the series increase, the methods for finding multiple change 
points do not perform well. Therefore we follow an iterative method for detecting change 
points in our special data. At the first stage, we apply a method for finding a single change 
point on our series. We use the non parametric method based on the cumulative sums 
test statistic (see [15] for details). The method returns the single most probable among all 
possible changepoint locations. Then we delete the data from this point to the end of the 



(3.4) as 



i+d* 




(4.1) 



6 



data, and again we iterate the method in the remaining data. 



Second Method: 

This method is based on calculation of sample variances of WiS, defined by ()4.ip . for i = 
1,2 - ■■ z and for i = z + 1, ■ ■ ■ n* as sample variance of those WiS preceding some point z as 
L{z) and those succeeding point z by U{z), as 

Liz) = - - Wi,,f, U{z) = - ^2,.)', 

z ^-^ n* — z ^-^ 

i=l i=z+l 
1 ^ - 1 "* 

= - E = - — 

z ^-^ n* — z ^-^ 

i=l i=z+l 

for all z = /*, /* + 1, . . . , n* — ^* where n* defined by (j4.1|) . and for some appropriate /* G N, 
say /* = 30. Then we evaluate 

S{z) = L{z) + U{z) for ah z = I* ,1* + 1, . . . ,7i* - I* (4.2) 

and plot them with respect to time z. As it is shown in Figure 4, the minimum of this plot 
occurs at the starting point of the last scale interval, say ii. By omitting samples from the 
point ii — j* to the end, and by repeating this technique for samples Wi, W2, ■ ■ ■ Wi-^-j*, 
for some appropriate j* € N, say j* = 50 one could find the minimum of S{k),S{k + 
1), . . . S{ii — j* — I*) at the starting point of the last scale interval of the remaining samples, 
say i2- Again by omitting samples from the point 12 — j* to the end and by repeating this 
method for the samples Wi, W2, ■ ■ ■ Wi^-j*, the starting of the last scale intervals appears 
as the minima of the remaining samples, as ^3. Then we evaluate the initial estimation of 
scale parameter, say Aq as ^E^- 

4.2 Estimation of scale parameter 

Assume that {X{t),t > 1} is a discrete scale invariance (DSI), also called semi-selfsimilar, 
process with some unknown scale A. We present an iterative method for the estimation 
of the scale parameter A and Hurst parameter H for the case that the process has sta- 
tionary increments. Successive Scale intervals for such process is considered as [A'^~^,A'^), 
where k € N. We also assume that the process could have self similar property with some 
prescribed Hurst index H' for samples inside each scale interval. Proposing equally space 
sampling implies that the increments have the same distribution. The iterativel estimation 
method is described for sfBm described in section 2.1 by the following steps. 

1- The whole duration of the study of the process is considered as the time interval [1,C]. 

2- Sampling of the process is assumed at points tQ,ti, . . . ,tN, where 

C-1 

ti — ti-i + 



for i = 1,2, . . . , N , and to = 1- 
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3- For implementing this method, we needs some initial value for the scale parameter, say 
Ao > 1, which is evaluated by a numerical method, described in section 4.1. 



4- Successive Scale intervals are considered as [Aq ^, Aq), A: = 1, . . . , M — 1, where M, the 
number of scale intervals, is evaluated as the largest / E N where Ag < C. 

5- Increments of the process are denoted by Y[ti) = X{ti) — for i = 1, . . . . As 
E[B}ji{ti)] = 0, so sample variance of the increments in A;-th scale interval for this initial 
scale value Aq, /c = 0, 1, . . . , M — 1 can be written as , S'^{\q) can be written as 



1 \2k{H-H') ^fc 

where = Yl!i=i "^ii — 0) ^^'^ ^« the number of samples in z-th scale interval. 

6- For this iterative method, several other points, say more than 50 points, are considered 
in each side of Aq and at some small equally distances of each other. Assume that the total 
number of such points round initial estimate Aq is denoted by m. Rename these ordered 
points as ai,a2, • • • flm- For all these points steps 4 and 5 are followed. Then for each a, 
we evaluate sum of the sample variances corresponding to some of the last scale intervals 
evaluated in step 5, say j number, which covers at leat %95 of observations as 

M 

R{a,)= Slia,). (4.3) 

k=M-j 

We plot such R{ai) , . . . , R{am) , which provide two clusters, those i?(ai)'s corresponding 
to the Oj's less than true scale parameter, and those corresponding to a^'s greater than 
true scale parameter. In example 1, such clusters are clearly shown by Figures 5(b) and 
5(e). Then for detecting change point of such identified clusters and estimate true scale, we 
evaluate sample variance of those R{ai) preceding each prescribed points as L*{ak) and 
sample variance corresponding to those R{\i) succeeding such prescribed as U*{ak), for 
all afc*, . . . ,a„i-k*j for some appropriate k* , say k* = 20. Thus 

L*{ak) = ^—^ J^(i?(ai) - Rif, U*{ak) = Yl (^(«^) " ^2)'- 

1=1 i=k+l 

where k = k* , . . . , m—k*, in which ^1 and R2 are corresponding sample mean of i?(ai), . . . , R{ak) 
and of i?(a/c+i), . . . , R{am) respectively. Then we evaluate and plot 

V{ak) = L*{ak) + U*{ak) (4.4) 

for these points. The point where V{-) has a minimum, which is considered as proper esti- 
mation of scale parameter is called A*. In example 1, plots 5(c) and 5(f) show the minimum 
of such V{-) for case 1 and case 2. 

7- Evaluate ftk = -J^^tttt, for A; = M - j, . . . , M. 
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8- Use the equation fi* = (A*)^^^ ^ \ to estimate H — H', where fx* is the weighted average 
of Hk] M — j<k<M — 1 evaluated at step 7, as 

AI M 
k=M-j k=M-j 

As one could estimate A* by step 6, this relation can be used to estimate H — H' . So it 
remains to estimate H' , the Hurst index of underlying fBm. By dividing samples of the 
A;-th scale interval to ) fQj- ]^ = 1,2,... successive observations of fBm can be 

evaluated. Even though one can follow the usual methods for estimation the Hurst index 
of self-similar processes, we present a new method for such estimation by the followings. 



4.3 Estimation of Hurst parameter for H-sssi processes 

In this section we present two new methods for estimating Hurst parameter of self-similar 
process with stationary increment (H-sssi), which can be followed to estimate Hurst index 
H' of the sfBm. Let {Xi,i = 1,2,..., N} be equally spaced samples of some H-sssi process 
as the main samples. We consider some sub-samples at points {Xj.fc, i = 1,2,..., [A^/A:]} as 
the k-th sub-sample for some fixed /c € N. The choice of the values of k depends on the 
sample size we take for instance k € {1, A'max}- 

For every k € {1, ATmax} we consider two sub-samples as {Xi} and {Xj fc}, where 
i = 1,2, ... , [N/k]. For these sub-samples, first and second order increments are defined 
as Yi^i = Xi+i - Xi and ¥2^1 = Aj+2 - 2Aj+i + Xi and as Yi^j.^ = A(j+i).fc - Xi.j, and 
^2,i fc = Ar(j_|_2).fc — 2X(j+i).fc respectively. By the followings we introduce two different 
methods for estimation Hurst index, which are designed to evaluate k"^^ by the ratio of the 
sample variances of the first and second order increments of such sub-samples respectively. 
The first method is accurate for H < 0.75 and the second method is more accurate for 
H > 0.75. 

First Method: 

For every k € {1, i^max}) we evaluate sample variances of first order increments of the 
above mentioned sub-samples. So we calculate sample variance of the first order increments 
as Yi^i, Yi^2; • • • ^i,[Ar/fc]-i5 as Sff^^. Also sample variance of the first order increments 
^i,fc)^i,2-fc5 ■ ■ ■ ) ^i,([Ar/fc]-i)fc as Sf 2) and evaluate H'k by the relation 

- e^'^ (4.5) 
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where 



q2 

'-'l,fc,2 — rAT 



UJ ^ i=i IfcJ ^ i=i 



and 



Yi,k,2 = tJ— E Yi,.k = j£— E ^1.^ = ^""'^M^i 

k\ ^ i=l LfeJ ^ i=l 
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Figure 3: Left three plots are 1-sfBm, 2-Increments of sfBm, and 3-Moving average Wi of MSV of the increments 
of sfBm, all with scale 2. The right three figure are corresponding figures with scale 4.. 



By ()4.5p we find H'l; and finally estimate H' as the mean of different H'/^ which have been 
evaluated as: 



ax 



k=2 

Second Method: 

We follow the same steps as mentioned in the first method, but this time based on the 
second order increments of sub-samples. So we evaluate sample variances of the second 
order increments 12,1, ^2,2, • • • ^2,[7V/fc]-2> ^-s S"! k i' ^^"^ sample variance of the second order 
increments l2,fc5 ^2,2 /0 • • • ) ^2,{[Ar/fc]-2)A: a-s 2^ evaluate H'k by the fact that 



02 

'^2,fc,l 



(4.6) 



where 



S. 



2,fc,2 rAT-l o 

^TJ i = l 



d k' 



2H' 



m-2 



in which 



2,k,2 



riVi 



J] {Y2,^ - Y2,kX 



2 _ i,2H' c2 
— '-'2,A;,lJ 



By (j4.6p we find W \^ and finally estimate H' as the mean of different ll\ which have been 
evaluated for different k by this method as: 



Example 1 

In this example we simulate sfBm processes defined by (2.3) with scales A and Hurst param- 
eters {H, H') as semi-selfsimilar processes with stationary increments, and then apply our 
methods to estimate their parameters. For doing this we consider two cases first, Ai = 2, 
Hi = 0.9, H'l = 0.2 and second, A2 = 4, = 0.6, H'2 = 0.2 and simulate 100,000 samples 
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Figure 4: 4(c)-Piot of Sum of sample variances 5(2), defined by 114.21 1 for the whole samples Wi, i = 1, - ■ ■ ,n*, for 
case 1, example 1, which has minima at the start of last scale interval, say ii. 4(b)- Plot of S{z) but for remaining 
samples Wi, when samples from 13 — I* to the end of data is removed, which has minima at 12, starting point of the 
last scale interval of remaining samples. 4(a)- Plot of S(z), this time by removing samples from 12 — j* to the end of 
data, which has minima at ii, starting point of the last scale interval of remaining samples. The last three figures are 
corresponding figures for case 2 of example 1. 



in each case, which have been plotted in Figures 3 or plan is to estimate scales first. So 
we follow to estimate some initial value by method one and two, described in section 4.2. 
Using the first method we detect change points in the corresponding MSV plotted in Figure 
3 at the starting points of last three scale intervals as oi = 65480, bi = 32745, ci = 16354 
for the first case and as a2 = 65497, 62 = 16349, C2 = 4060 for the second case. So the 
corresponding initial value of the scale parameter would be Ai = 1.99 and A2 = 3.99 by the 
first method. 



By the second method we detect the starting points of the last three scale intervals for 
case one as 65519, 32743, 16351, which is shown by Figures 4(c), 4(b) and 4(a), which corre- 
sponding statistics S{z) is defined by relation (j4.2p . Thus initial values of scale parameter 
for case one is evaluated as Ai = 1.9995. For case 2 the starting points of the last three 
scale interval is estimated as 65498, 16345,4056 which has been shown by Figure 4(d), 4(e) 
and 4(f), as the minima of S{z) defined by (I4.2p . Thus initial value of scale parameter is 
evaluated as A2 = 3.9998. Then following the iterative iterative estimation method we eval- 
uate V{ak) which defined by (|4.4p for 1000 equally spaced points in the interval [1.95, 2.05], 
round initial estimate Aq = 1.9995, which has been plotted by 5(c), and has minima at 
A^ = 2 as the best estimate for scale parameter in case 1,. Also be evaluating V{ak) for 
2000 equally space points in the interval [3.9,4.1] round initial estimate Aq = 3.9998, which 
has been plotted by (5.f) and has minima at A2 = 4 as the best estimate for scale parameter 
in case 2. We evaluated sum of sample variances of five last scale intervals for case one 
and for three last scale intervals for case 2 with estimated scale parameters A^ = 2 and 
A2 = 4 respectively, so that our scale intervals cover at least %95 of samples. Then we eval- 
uated fll = 2.6480 and — 3.0108 corresponding to \\ and A| respectively as described 
in step 8 of sequential estimation in section 4.1. Therefore we evaluate our estimations as 
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Figure 5: 5(a)-Cumulative sum of moving sample variances of increments of sfBm of case 1, examplel, 5(b)- 
Sum of sample variances of increments of R{ai), defined by 1 14.31 1. for 1000 equally spaced samples in the interval 
[1.9501, ■ ■ ■ 2.05] round the initial estimate of scale Ao = 1.9995, and 5(c)- Sum of sample variances V(ai^), defined by 
1 14.411 of preceding and succeeding samples to each point aj, in Figure 4(b), which has minima at the change point of 
figure 4(b), A* = 2, as the best estimate of scale parameter, for case 1, examplel. Right three figures are corresponding 
figures for case 2, example 1, that shows detection of the best estimator for scale as A* = 4 in Figure 5(f). 

Hi - H'l = 0.7024 and H2 - H'2 = 0.3992. 

Finally we extract the corresponding samples of underlying fBm as H-sssi process, by 
dividing samples in /c-th scale interval of sfBm in case 1 by {\\)^^~^ ^ = Ji\ = 1.6273, so 
Hi - H'l = 0.7025, and in case 2 by (A^)^^"^'^ = A2 = 1-7352, so H2 - H'2 = 0.3976. 
By applying our method described above in section 4.3, we find for case 1 that H'l = 
0.1957, Hi = 0.8981 and for case 2. So we evaluate H'2 = 0.2039 and H2 = 0.6031. 

Efficiency of the Hurst parameter Estimation 

To visualize efficiency of our Hurst estimation methods, we have simulated 10000 samples 
of fractional Brownian Motion with different Hurst parameter with 500 repetition and pre- 
sented the graph of the mean square errors (MSE) for Hurst Estimation of our methods and 
also of the method of Quadratic Variation [18] and the method of Convex Rearrangement 
|24j . We have used the difference of order two for the last two methods. As it is shown 
by Figure 6 for H < 0.75 our first method has less mean square error in compare to other 
methods and so is the best. For H > 0.75 our second method has less mean square error in 
compare with other methods and so is the best. 

5 Conclusion 

This paper provides some innovative method for estimation of the scale parameter of some 
class of semi-selfsimilar processes. For this we consider some starting point for iterative 
estimation of scale parameter and then by some precise method we estimate it. There are 
various open issues and vivid discussions about the estimation of scale parameter but still 
there is not a universal method which could be considered as the most promising method to 
find the best approximation of scale parameter in all cases. For the processes with station- 
ary increments, considering equally space sampling implies that the increments have the 
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Figure 6: Mean square error in Estimation of Hurst index of using 10000 samples of fBm with 500 repetition. 

same distribution inside each scale interval, so this enables us to consider such an iterative 
estimation method based on moving sample variances. Simulations and numerical evalu- 
ations clarified our results for the simple fractional Brownian motion. By detecting start 
points of consecutive scale intervals, we obtained an initial value for the scale parameter. 
Moving sample variances provide some clusters where jumps occurs at the starting point of 
scale intervals. Our iterative method provided a very good estimation for scale parameter. 
We also presented a heuristic method for the estimation of Hurst parameter of self-similar 
process which has less MSE in compare to previous methods. Moreover our estimation 
methods are easily implemented and are computationally fast. 
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